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ABSTRACT 

We have analyzed the data on 16,836 RR Lyrae (RR Lyr) variables observed toward the Galactic 
bulge during the third phase of the Optical Gravitational Lensing Experiment (OGLE-IH), which 
took place in 2001-2009. Using these standard candles, we show that the ratio of total to selective 
extinction toward the bulge is given by Rj = Ai/E{V — I) = 1.080 ± 0.007 and is independent 
of color. We demonstrate that the bulge RR Lyr stars form a metal-uniform population, slightly 
elongated in its inner part. The photometrically derived metallicity distribution is sharply peaked at 
[Fe/H] = -1.02 ± 0.18, with a dispersion of 0.25 dex. In the inner regions {\l\ < 3°, \b\ < 4°) the RR 
Lyr tend to follow the barred distribution of the bulge red clump giants. The distance to the Milky 
Way center inferred from the bulge RR Lyr is Rq = 8.54 ± 0.42 kpc. We report a break in the mean 
density distribution at a distance of ~ 0.5 kpc from the center indicating its likely flattening. Using 
the OGLE-HI data, we assess that (4-7) x 10^ type ab RR Lyr variables should be detected toward 
the bulge area of the on going near-IR VISTA Variables in the Via Lactea (VVV) survey, where the 
uncertainty partially results from the unknown RR Lyr spatial density distribution within 0.2 kpc 
from the Galactic center. 
Subject headings: Galaxy: bulge - Galaxy: structure - stars: variables: other (RR Lyrae) 



1. INTRODUCTION 

RR Lyrae (RR Lyr) variable stars are known as stan- 
dard candles and useful tracers of old populations. They 
play an essential role in our understanding of the for- 
mation and evolution of the Galaxy. When multi-epoch 
light curves are available, RR Lyr variables can be easily 
identified. A large number of observations allows deter- 
mination of their properties with high precision. The 
derived parameters can be used to determine interstellar 
extinction, metallicity, distance, and spatial distribution 
of the stars. RR Lyr stars observed toward the bulge 
provide an independent determination of the distance to 
the center of the Milky Way. 

RR Lyr stars observed close to the central regions of 
the Milky Way co ncentrate toward the Galactic center 
was first noticed bv lvan GentI ()1932l[l"933l ). Later. iBaadd 
([1946) found a strong predominance of RR Lyr vari- 
ables toward a relatively unobscured area today called 
the Baade's window, indicating the presence of Popula- 
tion II stars in the Milky Way center. By the early 1990s 
approximately 1000 RR Lyr variables inhomogeneously 
distributed toward the bulge were known. 

The number of new bulge RR Lyr variables has in- 
crease d following the advent of mas siv e pho t ometr ic sur- 
veys. lUdalski et all (IT991 [T995al lhl. 119961 [l99l pub- 
lished a catalog of over 3000 periodic variable stars de- 
tected in the fields covered by the first phase of the Op- 
tical Gravitational Lensing Experiment (OGLE-I). Two- 
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hundred and fifteen of these objects were classified as 
RR Lyr. Analysis of the O GLE-II d ata brought a much 
larger list of 2713 RR Lyr (Mizerski' 200 ^). La ter, using 
the same source of data, iCollinge et al.l (|2006[ ) prepared 
a catalog of 1888 fundamental mode RR Lyr stars (type 
RRab). 

The MACHO microlensing project also observed a nu- 
merous sample of RR Lyr stars toward the Galactic 
center. iMinniti et al.l (|1998| ). using a sample of 1150 
RRab and 550 RRc stars, showed that the spatial dis- 
tribution of RR Lyr between 0.3 and 3 kpc follows a 
power law with an inclination of —3.0. lAlcock et al.l 
(1998) examined the mean colors and magnitudes of 
^1800 RR Lyr and found that the bulk of the pop- 
ulation is not barred. Only RR Lyr located toward 
the inner fields closer to the Galactic center {I < 4°, 
b > —4°) seem to follow the ba rred distribution ob served 
for red clump giants (RCGs) fStanek et al.l ll994[ ). Re- 
cently, Kundcr ct al. (2008) analyzed photometric data 
on 3525 MACHO RRab stars to assess the reddening to- 
ward the Galactic bulge. They derived the selective ex- 
tinction coefiicient Rv,vr = Av/E{V - i?) = 4.3 ± 0.2, 
which corresponds to the average value observed in the 
solar neighborhood Ry.BV = Av/E(B -V) = 3.1 ± 0.3 
(|Cardelli et al . I [l989tlFiFzpatrickl 119991) . 

The first metallici ty measurem e nts of bulge RR Lyr 
stars wer e made bv iButler et all ()1976[ ) using the AS* 
method (jPrestonl I1959D. Us"ing 9 sta rs in the 100 
variable sample of iBaade et al.l (|1963D list, they ob- 
tained ([Fe/H]) — —0.65 ± 0.15 dex and concluded 
that t he stars are mildly metal-poor. iGratton et al.l 
(jl986f ) used a sample of 17 bulge RR Lyr variables and 
found a wide range of iron abundances, between —1.8 
and +0.1 dex. Later, from spectr a of 5 9 RRab and 
RRc variables, iWalker fc Terndrunl (|199lD determined 
an average metallicity of {[Fe/H]) = —1.0 dex on the 
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iZinn fc Westi ()1984D metallicity scale. In contrast to pre- 
vious studies they found the metallicity distribution to 
be sharply peaked, with a dispersion of only 0.16 dex. 
The most recent metallicit y estimate was performed by 
iKunde r fc Chaboved ()2008f ). Based on photometric data 
for 2,435 MAC HO RRab stars th ey derived ([Fe/H]) 
= -1.25 (on the iZinn fc Westlll984l scale), with a broad 
metallicity range from [Fe/H] = —2.26 to —0.15 dex. 

The bulge RR Lyr stars have also often been used 
to measure the distance to the Galactic center, un- 
der the assumption that the center of their popu- 
lation corresponds to the cente r of the Mil ky Way. 
The first such measurement by iBaadd (|1951[ ) yielded 
Rq — 8.7 kpc. For more than half a century, vari- 
ous investigations brought d ifferent value s of Rq with 
comparable uncertainty, e.g. . IHartwick et al. (1972:) ob- 
tained Rq = 7. ± 0.6 k pc, Oort & PlautI (|1975l ) Rq = 
8.7 ± 0.6 kpc, iBlanco & Blanco (1985) Rq = 6.94 ± 
0.58 kpc, [W alker & Mack (1986) Rq = 8.1 ± 0.4 kpc, 
iFernlev et al., (,1987, ) Rq = 8.0±0j5 kgc^ lCarnev etaLl 
( 19951) R q = 7.8 lb 0.4 k pc. iFeast I (119971) Rq = 8.1 ± 
0.4 kp c, CoUingc ct al.jj2006D Rq = 8.3 ± 0.7 kpc, 
iGroencwcgcn ct al. (20 0l Rq = 7.94 ± 0.63 kpc, and 
recently [Mijaesd (|2010l ) i?o = 8.1 ± 0.6 kpc. 

In this paper, usi ng the newly released catalog 
(jSoszyhski et all l2011f) of nearly 17,000 RR Lyr stars 
detected toward the Galactic bulge during the third 
phase of the OGLE project (OGLE-HI, Udalski 2003a; 
lUdalski et al.l l2008') , we show that these variables consti- 
tute a uniform population, different to other populations 
residing in the central regions of our Galaxy. 

The structure of this paper is as follows. Section [2] de- 
scribes the selection of the sample and gives information 
on the reddening. Metallicities and magnitudes of the 
bulge RR Lyr variables are determined in Section [3l In 
Section m we estimate the distance to the Galactic cen- 
ter. The structure of the RR Lyr population, based on 
their mean brightness and density distribution, is studied 
in Section [5l Section |6] states our main conclusions. 

2. THE DATA 
2.1. Cleaning the Sample 

iSoszvnski et al.l ()2011[ ) published data on 11,756 RRab, 
4989 RRc, and 91 RRd stars observed toward the Milky 
Way bulge. The mean /-band magnitudes of the detected 
stars lie in the range 12 < / < 20. For the purpose of 
our work, the sample had to be cleaned from various 
contaminants. In the first step we rejected RR Lyr with 
amplitudes lower than 0.08 mag in /, which are very 
likely blended objects. Figure [1] presents /-band ampli- 
tude distributions for each of the three RR Lyr variable 
types. Based on the shape of the distributions for RRab 
and RRc stars we set the lower limit at an amplitude 
of 0.08 mag resulting in the rejection of 38 RRab, 17 
RRc, and RRd objects. We also removed the object 
OGLE-BLG-RRLYR-02792, which likely belongs to a bi- 
nary system, and a RRab variable with a suspiciously 
large amplitude of 2.0 mag. The remaining amplitudes, 
of stars left in the sample, were no higher than 1.49 mag. 

In the next step, we constructed / versus V — I dia- 
grams to separate bulge RR Lyr stars from foreground 
and background objects. The diagrams for RRab, RRc, 
and RRd variables are shown in FiguresIlHll respectively. 
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Fig. 1. — Amplitude distributions for all three types of RR Lyr 
varia bles detected in the OGLE-III bulge area (Soszvhski et al] 
120111) . A bin size of 0.02 mag is used. Stars with 7-band am- 
plitudes lower than 0.08 mag (marked with a thin vertical line) are 
very likely blended objects and were rejected from the analysis. 
Note that the histograms show logarithm of the counts, thus bins 
with a single star appear blank. 

In all these figures one can see the effect of heavy redden- 
ing. Most of the background RR Lyr stars belong to the 
Sagittarius dwarf spheroidal (Sgr dSph) galaxy. This is 
well seen in the case of numerous and bright RRab stars 
which form a sequence in the lower part of Figure [2] 

Finally, from the list of bulge RR Lyr stars we rejected 
51 RRab and 23 RRc variables which are members or 
likely members of eight Galactic globular clusters (GCs; 
P. Pietrukowicz et al., in preparation). 

2.2. Reddening 

The selection returned 10,472 RRab, 4608 RRc, and 78 
RRd bulge variables with measured mean V- and /-band 
magnitudes^ The corrected /-band magnitudes are ac- 
curate to 0.01-0.02 mag for {V — I) < 6 mag. For an 
additional 627 RRab, 151 RRc, and 2 RRd stars that 
are located close to the Galactic plane, we do not have 
information on l^-band brightness due to the heavy red- 
dening in those directions. We used the information on 
the /-band brightness and V — I color to estimate the 
ratio of total to selective extinction Rj = Ai/E{V — I), 
and to test whether it is independent of color or not. In 
Figures [S] and [Bl respectively for RRab and RRc stars, 
we show unbinned and color-binned bulge RR Lyr stars 
in the / versus V — I diagram. The observed scatter 
of ^ 0.4 mag is due to intrinsic differences in the abso- 
lute brightness, differences in extinction along the line of 
sight to the RR Lyr star, the effect of the extent of the 

^ We note t h at al l magnitudes published in the ca talog of 
'Soszyiiski ct al. ( 20111 ) were calibrated to the standard Landolj 
(1992) photomet ric system using transformation formulae pre- 
sented in lUdalski et al .l (2002) a nd an improved relation for stars 
with (V - I) > 2 mag in Szvmahski et al.f (|20111 'l. 
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Fig. 2. — Color-magnitude diagram (based on mean values) for 
11,124 OGLE-III RRab variables. For further analysis only stars 
located between the two inclined lines were used. Most of the 
background objects belong to the Sagittarius dwarf spheroidal (Sgr 
dSph) galaxy. Accuracy of the mean brightness is estimated as no 
worse than 0.02 and 0.05 mag in the / and V bands, respectively. 
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Fig. 3. — Same as in Figure [2] but for 4837 RRc variables. 



population along the line of sight, and the photometric 
uncertainty. We mark the limits inV — I between which 
we fit a linear regression. Stars with high V — I are very 
reddened and too faint to be observed. On the other side, 
there are very few mildly reddened bulge RR Lyr stars. 
We obtained the following relations for RRab and RRc 



Fig. 4. — Same as in Figures [2] and [3] but for 89 RRd variables. 

stars, respectively: 

/RRab = (1.085±0.012)(y-J)+(14.443±0.021), a = 0.42, 

(1) 
/rr^ = (1.075±0.014)(y-/)+(14.712±0.022), a = 0.36. 

(2) 
From the fits we find a mean value of Rj — 1.080±0.007. 
We also find the difference in mean /-band brightness be- 
tween RRc and RRab stars to be A/ = 0.27 ± 0.03 mag. 
The obtained value of the ratio of total to selective ex- 
tinction is in agreemen t with Rr = 1. 1 inferred from 
the OGLE-II data by lUdalskil fl2003bn FI Our result 
indicates that interstellar extinction toward the bulge 
is indeed anomalous if compared with the standard 
value of Ri.vi = 1. 4 derived from the all-sky maps of 
ISchleeelet aLllfTOOl) . 

2.3. Selection of the Subfields 

The field of view of the OGLE-III 8-chip camera was 
a square 35f6 on a side. We used this area as a unit 
field in our analysis. To obtain good statistics of RR 
Lyr stars, we divided the unit field into four and fur- 
ther into sixteen square subfields, corresponding to two 
adjacent OGLE-III chips (17f8 on a side) and half of a 
chip (8.'9 on a side), respectively. The lowest number of 
RR Lyr stars required per subfield we set at 10. If the 
number of stars was lower than 10, we counted stars in 
a larger subfield. Poorly populated unit fields (with less 
than 10 stars) were not taken into account in the analy- 
sis of mctallicity and brightness distributions (Section[3l). 
In Figure [71 we show a map of the bulge RRab stars in 
Galactic coordinates together with marked centers of se- 
lected subfields. The number of RRab subfields is 396 



6 A later determined value of 1.168 bylNataf et al.l (201(1) from 
brightness of OGLE-III red clump stars is slightly different, since 
they calibrated photometry for stars redder than {V — I) = 2 mag 
using an older and slightly inaccurate formula. 
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Fig. 5. — Color-magnitude diagram with a linear fit for unbinned 
(upper panel) and binned (lower panel) bulge RRab variables. The 
bin size is 0.1 mag. Thin vertical lines mark section to which we 
fit the straight line. 



Fig. 7. — Location of 11,099 bulge RRab variables in Galactic 
coordinates. Large points denote centers of 396 square subfields. 
Horizontal lines limit 10 selected latitudinal stripes used in further 
analysis. With the rectangle we ma rk the bulge area of the VVV 
near-IR survey ((Minniti et an i2010l '). 
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Fig. 6. — Same as in Figure [6] but for RRc type variables. Note 
very similar value of Rj for RRab and RRc stars. 

with 10-64 stars per subfield. In case of the RRc stars 
we selected 200 subfields, each containing 10-53 stars. 

2.4. Completeness of the Sample 

Completeness of the OGLE-III catalog of RR Lyr stars 
depends on many factors such as brightness of stars, their 
amplitude, shape of the light curve, crowding, and num- 
ber of observing points. The search for periodic variabil- 



ity was run on /-band light curves with more than 30 
points for about 3 x 10* stars with brightness down to 
I ^ 21 mag. The light curves typically consist of 100- 
3000 measurements. In order to test and impro ve com- 
pleteness of the catalog, iSoszvhski et al.l (|2011[ ) carried 
out a cross-identification with publis hed catalogs o f RR 
Lyr stars base d on MACHO (Kund er et "all 120081 ) and 
OGLE-II data (|Collinge et al..,2006, ). Out of 2114 MA- 
CHO RRab stars they missed only 27 objects. Thirteen 
of them turned out to be located close to bright, sat- 
urated stars that were masked during reductions. The 
cross-match with 1888 OGLE-II RRab variables resulted 
in one missing object, a blended star. Based on the cross- 
identification, iSoszvnski et al.l (|2011[ ) estimated the com- 
pleteness of the bulge RRab variables at a level of about 
99%. 

In Figure [8l we compare the luminosity histograms 
for all detected OGLE-III RRab and RRc variables 
with the histograms for all stars detected stars toward 
three different Galactic latitude intervals. The OGLE 
photometry is complete down to / ^ 19.5 mag and 
there are very few RR Lyr stars with mean magnitudes 
fainter than this brightness. This shows that the limit- 
ing magnitude weakly affects the total number of de- 
tected RR Lyr stars. The completeness, at least for 
bright RRab stars with characteristic high-amplitude 
saw-tooth-shaped light curves, is indeed very high. RRc 
stars, on the other hand, are on average fainter than 
RRab variables and have low-amplitude sinusoidal light 
curves that are very similar to the light curves of ellip- 
soidal binaries. Hence, pulsating stars of this type are 
more difficult to classify. 

3. METALLICITY AND BRIGHTNESS OF THE BULGE RR 
LYRAE STARS 
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Fig. 8. — Luminosity functions of bulge RRab (dotted line), RRc 
(thin solid line), and all stars (thick solid line) detected in the 
OGLE-III /-band data within 1?0 around {l,b) = (0,-2) (upper 
panel), {l,b) = (0,-4) (middle panel), and {l,b) = (0,-6) (lower 
panel). Note that the OGLE photometry is complete down to 
/ ~ 19.5 mag, which is deeper than the mean brightness of almost 
all detected RR Lyr stars. 



3.1. Metallicity 

According to the idea of iKovacs fc Zsoldo^ (|1995| ) , the 
iron abundance of RRab stars can be derived from their 
light cu rve parameters. Th e method was further devel- 
oped bv lJurcsik fc Kovlcsl (|l996) who presented a basic 
relation between the metallicity [Fc/H], pulsational pe- 
riod P, and Fourier phase combination c/f^i — (j)3 — 30i 
from V^-band light curves. We determine metallicities 
for 10,259 reliable bulge RRab variables using a relation 
suitable for J-band light curves found bv lSmole"3 ()2005l ) 
for sine decomposition: 

[Fe/H]jg5 = -3.142 - 4.902P + 0.824031, asys = 0.18. 

(3) 
The derived metallicities are on the iJurcsikl (jl995^ scale 
(J95), which is based on high dispersion spectroscopic 
measurements. We used the uncertainty of (f>3i to cal- 
culate the intrinsic metallicity uncertainty a-mt of each 
star. They range from 0.002 to 0.30 dex with a median 
value of 0.012 dex. This uncertainty quadratically added 
to cTgys given by Smolcc (2005) results in a median total 
uncertainty of crpc/H] — 0.18 dex. 

In Figure [9] we show the metallicity distribution for 
the bulge RRab stars, plotted as a histogram with a bin 
size of 0.05 dex. In the upper panel of this figure we 
plot, for comparison, a distribution of average RR Lyr 
metallicity for the 396 subfields which were analyzed. 
Both distributions are sharply peaked, which indicates 
that the bulge RR Lyr population was likely formed 
on a short timescale, as it is in the case of t he Galac- 
tic halo and thick disk ([Gratton et al.l I1996D . About 
98.7% of the bulge RRab stars have metallicities be- 
tween —2.0 and 0.0 dex, with a peak at [Fe/Hjjgg = 



— 1.02 and a dispersion of 0.25 dex. These numbers 
are in very good agreeme nt with results obtained by 
iKunder fc: Chaboveij (|2008l) from F-band light curves of 
2435 MACHO RRab variables. They found an aver- 
age metallicity of ([Fe/H]) = —1.25 and a dispersion 
of 0.30 dex on the Zin n fc West (1984) metallici t y scal e 
which, according to the work of iPapadakis et all (|2000D , 
is less metal- rich by about 0.24 dex than the scale of 
IJurcsikl (Il995l) . 

In the lower panel of Figure |9l we compare the bulge 
RRab mctallicity-distribution function derived in this 
work with that of the 47 Galactic globular clusters (GCs) 
located within 3.0 k pc from the G alactic centeiQ (based 
on 2010 version of iHarrisI 119961 catalog) and that of 
a sample of 403 bulge red giants (RGs) published by 
iZoccali et all ([200.80. The metallicities of the GCs and 
RGs were converted from th e lCarretta et a l. (2009) scale 
(C09) to the one of IJurcsikl (|1995[ ) via the .Zinn fc West! 
(I1984D metallicity scal e (ZW84) using a relation from 
IPapadakis et al.l|2000( ): 

[Fe/H]zw84 = l-028[Fe/H]j95 - 0.242, (4) 

and a second relation from iCarretta et al.l (|2009| ): 



Fe/H]co9 = 1.105[Fe/H] 



ZW84 



0.160. 



(5) 



The metallicity distributions of the bulge RGs, RR Lyr 
stars, and GCs differ significantly from each other. The 
very peaked and symmetric distribution for the RRab 
stars indicates that the population of the bulge RR Lyr 
variables is very metal uniform. In Figure [TOj we show 
a color-coded map of metallicity averaged in each of the 
396 subfields. The metallicity of RR Lyr stars seems to 
be fairly uniformly distributed over the whole investi- 
gated area. 

3.2. Dereddened Mean Brightness 

We calculate dereddened mean magnitudes /q for each 
RR Lyr star as: 

/q = / - Ai, 



(6) 
(7) 
(8) 



where the extinction is given by 

Ai=RiEiV-I), 
with reddening given by 

E{V-I) = {V^I)-{V-I)o. 

We compute the intrinsic color {V — I)o for RRab stars 
as 

{V - /)o = Mv - Mi, (9) 

where the absolute bright nesses My and Mr c ome from 
theoretical calibrations in iCatelan et al.l (j2004f ) 

My == 2.288 -I- 0.882 log Z + 0.108 (log Z)^ (10) 

M/ = 0.471 - 1.132 logP + 0.205 logZ, (11) 



^ According to IClement et al.l 11200 II ) and the updated Cata- 
logue of Variable Stars in Galactic Globular Clusters (2011) by 
C. Clement (http://www.astro.utoronto.ca/~ccloment/read.html) 
the 47 bulge (JUs contain d'Zb known RK, Lyr stars considered to 
be members of the clusters. Around two-thirds of those variables 
(exactly 217 stars) belong to the globular cluster M62. There are 
no known RR Lyr stars in 31 of the selected bulge GCs. 
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Fm. 9. — Normalized metallicity distribution (on the IJurcsikl 
11199 5') scale) for 10,259 bulge RRab stars (solid thick line) in com- 
parison with the distribution of average RR Lyr metallicity for the 
396 subfields that were analyzed (upper panel) and in comparison 
with distributions for 47 bulge GCs located within 3.0 kpc from the 
Galac tic center, and 403 bulge red giants (RGs) from the catalog 
of Zo ccali et al.l (p008) (lower panel). The distribution of RR, Lyr 
stars is sharply peaked at [Fe/H] = —1.02 dex and significantly 
differs from the populations of the bulge RGs and GCs. 
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Fig. 10. — Metallicity map for RRab variables averaged in 396 
analyzed subfields. The metallicity appears to be uniformly dis- 
tributed over the observed part of the Galactic bulge. 
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Fig. 11. — Distance distribution for 10,259 bulge RRab stars. 
The bin size is 0.2 kpc. The distribution is very well represented 
by a Lorentzian with a maximum value at Rg = 8.54 kpc. 

with the fohowing conversion between Z and [Fe/H] : 

logZ = [Fc/H] - 1.765. (12) 

The conversion is based on a solar metalhcity of Zq — 
0.01716, as required to match the m etal-to-hydrogen ra- 
tio Z/X of IGrevesse fc Noelsl (J1993I ). Intrinsic colors of 
RRc s tars are between .32 and 0.42 mag (based on mod- 
els in iFeuchtingeH Il999f ) . Since the metallicites of this 
type of RR Lyr variables are not determined, we use the 
mean value of {V — I)o = 0.37 mag for each RRc star. 

In the above calculation, we adopt an uncertainty of 
the mean /-band br ightness at a level of aj = 0.02 mag 
(|Udalski et al.l[2008h . The OGLE-III T/-band hght curves 
are less well sampled and the adopted accuracy of the 
mean I^-band brightness is try = 0.05 mag. The er- 
rors propagate to a range of ct/q — 0.07-0.20 mag for 
the sample of 10,259 RRab stars with a median value of 
0.09 mag. 

4. DISTANCE TO THE GALACTIC CENTER 

We estimate the position of the Galactic center by mea- 
suring the mean distance to the bulge RR Lyr variables. 
The distances R to individual bulge RRab star are esti- 
mated using the relation: 

logi?=l-F0.2(/o-M/). (13) 

To obtain the real distance to the center of the pop- 
ulation we have to apply the following two geometric 
corrections. First, we have to project all the individual 
distances onto the Galactic plane by taking the cosine of 
the Galactic latitude, yielding Rcosb. Second, we have 
to take into account the "cone-effect" — our subfields 
subtend solid angles on the sky with more volume fur- 
ther away — by scaling their distance distribution by 
R^^. A corrected and normalized distance distribution 
for the 10,259 RRab stars is presented in Figure [TT] The 
symmetry of this distribution further co nfirms the high 
compl eteness of the RR Lyr catalog of Soszvnski et al.l 
(|2011h . The slightly flat peak in the distribution between 
8.2 and 8.8 kpc is consistent with the presence of an elon- 
gated inner structure (see Section [5]). The distance dis- 
tribution is very well represented by a Lorentzian with 
a maximum value at 8.540 ± 0.013 kpc. A resulting dis- 
tance to the Milky Way center of Rq = 8.54 ±0.42 kpc is 
obtained by quadratically adding the median uncertainty 
in the distance to the individual stars of 0.418 kpc. This 
result is in good agreement with the weighted average of 
8.15 ± 0.14 ± 0.35 kpc determined from eleven different 
measuring methods bv .Gcnzcl ct al., (,2010.) . 
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Fig. 12. — Dereddened /-band magnitude distribution for bulge 
RRab stars (blue points), in comparison with red clump giants 
(RCGs, red points) in 10 latitudinal stripes (marked in Figure (TJl. 
Not shown are the individual errors on the magnitudes of RCGs, 
which are between 0.10 and 0.40 mag, with a median value of 
0.20 mag. The RR Lyr stars only follow the distribution of RCGs 
toward the inner part of the Galactic bulge (—3° < I < +3°, 
-4° < 6 < -2°). 



5. STRUCTURE OF THE BULGE RR LYRAE POPULATION 
5.1. Brightness Distribution 

In this section, we study the structure of the bulge RR 
Lyr population using their mean dereddened magnitudes 
at different Galactic latitudes. We selected 10 latitudi- 
nal stripes, they are 0?5 or 1?0 wide. Their location 
is illustrated in Figure [T] Figure [T^ shows Jq magni- 
tudes averaged in subfields for RRab stars in comparison 
to the mean brightnesses of RCGs from D. M. Nataf et 
al. (2012, in preparation). A typical uncertainty in the 
dereddened brightness of the RCGs is 0.20 mag. Results 
for RRc stars arc shown in Figure [T31 Bulge RR Lyr 
stars follow the barred RCGs distribution in the inner 
part of the bulge at |;| < 3° and |5| < 4°. This is weh 
seen in the better-covered southern part (see Figure [H)) . 
Farther off the Galactic plane at latitudes \b\ > 4° the 
distribution of RR Lyr variables is clearly flat, indicat- 
ing that their population becomes more spherical with 
increasing distance from the center. 

In Figure 1151 we investigate if there are any dif- 
ferences in dereddened magnitudes between metal-rich 
([Fe/H],g5 > -1.0) and metal-poor ([Fe/Hjjgg < -1.0) 
RRab stars. For each latitudinal stripe we fit a weighted 
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Fig. 13. — Same as in Figure [121 but for the bulge RRc stars. 

linear regression to the data. All derived mean values 
are positive, as expected from Equation ()lip . that metal- 
poor stars are on average brighter than metal-rich ones. 
All values are around -1-0.1 mag, indicating no significant 
differences in metallicity distribution along the line of 
sight of the bulge RR Lyr population for different Galac- 
tic latitudes. Moreover, there seems to be no significant 
trends in A/q with Galactic longitude. These results cor- 
roborate the metal-uniform nature of the bulge RR Lyr 
population. 

5.2. The Inner Part 

The OGLE-III V- and J-band data cover Galactic lat- 
itudes higher than b ^ -|-2?1 and lower than b ~ — 1?2. 
Using RRab variables with derived individual distances, 
we can make an attempt to study the shape and orien- 
tation of the observed part of the inner bulge. Figure [16] 
shows four projections of the variables onto the Galactic 
plane. Two of those projections, —3° < b < —4° and 
—6° < b < —5°, include complete data for stars with 
Galactic longitudes -4° < / < 4-4°. Figure [T7I illustrates 
the projection for —3° < 6 < —4° as a color-density map. 
The presence of the elongated structure for the RR Lyr 
located closer to the Galactic plane is evident. Farther 
off the plane the structure gets round, as is expected 
from the flat mean magnitude distribution presented in 
Figures [T2HI11 Based on the density map in Figure [T71 
we find the inclination of the RR Lyr structure major 
axis to be about 30° with respect to the line of sight 
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Fig. 14. — Comparison of the inclinations a for RRab stars (blue 
points) and RCGs (red points) in relation Ig vs. / for the Galactic 
longitudes —3° < Z < +3° in three well-observed latitudinal sec- 
tions. Uncertainties of the inclinations are given in parenthesis. In 
the inner regions (|fet < 4°) RR Lyr stars more closely follow the 
trend of the RCGs. 

between the sun and the Galactic center. This angle is 
very similar to the i nclination of 24°-27° der ived for the 
main RCGs bar bv IRattenburv et al.l ()2007[ ). However, 
we stress that the accurate estimation of the inclination 
of the RR Lyr spatial distribution is hampered due to 
small number statistics and incomplete coverage of the 
central part of the bulge. Near-IR data from the on go- 
ing VISTA Variables in the Via Lactea (VVV) survey 
(jMinniti et al. 2010) should bring the final answer to the 
question on the structure and shape of the inner bulge 
RR Lyr population. 

5.3. Search for Split in Brightness of the RR Lyr Stars 

Recently, iNataf et al.l (|20100 used OGLE-III photome- 
try to find that the Galactic bulge red clump is split into 
two components toward the fields at (—3? 5 < I < 1°, 
b < -5°) and (l,b) = (0 . +5?2 ). The split is also re- 
ported bv lGonzalez et ahl (1201 1[) who an alyzed the first 
data from the VVV survev (jSaito et al.| [2012'). The ob- 
served difference in /-band brightness between the com- 
ponents is approximately 0.5 mag. We selected OGLE 
fields located at |6| > 5° for which the number of detected 
bulge RRab stars was at least 20. Figure [18] presents 
luminosity histograms for the investigated fields. The 
adopted size of the bins, 0.25 mag, is wider than the 
typical brightness uncertainty of ct/q = 0.08 mag. The 
number of bulge RR Lyr variables is insufficient to con- 
clude if such a split is present in this type of stars or not. 
Some histograms display double peaks, but this can be 
explained by small number statistics. 

5.4. Density Distribution 

We analyzed the mean spatial distribution of the bulge 
RR Lyr as a function of distance from the Milky Way cen- 
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Fig. 15. — Differences in dereddened /-band magnitudes between 
metal-rich ([Fe/H] > -1.0) and metal-poor ([Fe/H] < -1.0) RRab 
stars in 10 latitudinal stripes. Mean differences (A/q) for all stripes 
are around a similar value of ~ -1-0.1 mag and there are no signif- 
icant trends in AIq with Galactic longitude, as is expected for a 
metal-uniform population. 



ter by counting the stars in square subfields that are 17f8 
on a side. Separate results for the whole sample (consist- 
ing of all RR Lyr pulsators), RRab, and RRc subsamples 
are presented in three panels of Figure[T21 The data cover 
angular distances between d « 1?1 and d « 11? 8 from the 
center, however, the innermost subfields are incomplete 
due to high reddening, while the outermost ones suffer 
from small number statistics. We found a break in both 
distributions at d w 3?5. For all detected RR Lyr vari- 
ables in the inner part between d = 1?5 and d = 3?5 we 
obtained the following fit: 



logs. 



allRR.inn 



(-1.13±0.08) logd+(3.28±0.04), (14) 



where d is in degrees and S in counts dcg^^. The outer 
part for 3? 5 < d < 8?0 is represented by a steeper line 
with 

logSaiiRR^out = (-1.66±0.06)logd+(3.57±0.04). (15) 
In the case of RRc stars we found, respectively, 

log SRRcinn = (-0.99 ± 0.13) logd+ (2.69 ± 0.06) (16) 
and: 

logSRRc,out = (-1.86±0.10)logd+(3.18±0.07). (17) 
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Fig. 16. — RRab variables from four different latitudinal stripes 
projected onto the Galactic plane, where the Sun is located at 
{x,y) = (0, 0). For stars located closer to the plane the elongated 
structure is evident. It disappears farther off the plane. Note that 
the coverage of the bulge area at some longitudes is incomplete (see 
Figure O 

Finally, for RRab stars we have 

logSRRab,i„„ = (-1.25±0.09)logd+(3.17±0.04) (18) 
and 

log ERRab.out - (-1.62±0.07)logd+(3.37±0.05). (19) 
This would correspond to a spherical distribution with 
i^RRab.inn oc r-2-25±o°9 for 0.2 < r < 0.5 kpc 
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-DRRab.out oc r 



-2.62±0.07 



for 0.5 < r < 1.2 kpc. 



assuming the distance to the Galactic center of Rq = 
8.54 kpc given in Section H) 

The distribution of RR Lyr stars in the central part 
(r < 0.2 kpc) remains unknown. We may speculate that 
it probably gets more flat toward the center. The an- 
swer may come from jjT.s-ban d observations curre ntly col- 
lected by the VVV survey (jMinniti et al.ll2010[ ). Using 
the OGLE-III data, we are able to assess the number of 
RR Lyr stars which should be detected within the VVV 
bulge area (marked in Figure [7]). If we extrapolate the 
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Fig. 17. — Density map of RRab stars with latitudes —4° < b < 
—3° projected onto the Galactic plane. Each bin is 0.12 kpc on a 
side. The inner structure is elongated and inclined by ~ 30° with 
respect to the line of sight. 
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Fig. 18. — Number histograms of RRab stars in selected subfields 
around Galactic latitudes -f 5° (upper row), —5° (upper middle and 
middle rows), and —6° (two bottom rows). Galactic coordinates 
for the centers of the subfields are given in each box. The size of 
the bins is 0.25 mag, i.e., almost three times wider than a typical 
photometric error (cr/o ^ 0.09 mag). Small number statistics pre- 
vent us from concluding whether a split such as that observed for 
RGG exists for RR Lyr. 
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Fig. 19. — Mean surface density distribution as a function of 
angular distance from the Milky Way center for all bulge RR Lyr 
(upper panel), RRab (middle panel), and RRc variables exclusively 
(lower panel). Three vertical lines at logd = 0.176, 0.544, and 
0.903, corresponding, respectively, to linear distance r = 0.2, 0.5, 
and 1.2 kpc, mark two sections of the data to which lines were fit. 
The outliers at logd « 0.25, 0.42, and 0.51 were not included in 
the fitting for the density relations. 



inner distribution toward the center, we find that VVV 
should detect (60 ± 7) x 10^ RRab stars. With a flat 
distribution in the central part we obtain (49 ± 7) x 10'^ 
stars. For RR Lyr variables of all types the numbers are 
the following: (82 ± 9) x lO^ and (71 ± 9) x 10^. In sum- 
mary, if we assume that the OGLE-III catalog of RRab 
variables is 99% complete, we can estimate the VVV sur- 
vey should detect (4-7) x 10^ RR Lyr stars of this type. 
In Figure [20] we draw density distributions of bulge 
RRab variables in longitudinal and latitudinal slices, 
each with width of 6°, to verify if there are any signifi- 



FlG. 20. — Mean surface density distributions for RRab stars 
of all metallicities (upper panel), metal-poor ([Fe/H]jgg < —1.0, 



middle panel), and metal-rich ([Fe/H]j 



> 



-1.0, lower panel) 



cant differences between metal-poor ([Fe/Hjjgg < —1.0) 
and metal-rich ([Fe/H]jgg > —1.0) stars. No differences 
are found in the linear fits, leading us to the earlier con- 
clusion that the RR Lyr bulge population is metal uni- 
form. The observed slope of the latitudinal distributions 
is always shallower t han the longitudinal one, confirming 
earlier suggestions of lMinniti et al.l (|1998D that the whole 
population is slightly flattened along Galactic longitude. 
Based on the longitudinal and latitudinal fits obtained we 
find a mean density of logE — 2.2 at distances diong = 4?9 
(logdiong = 0.69) and diat = 6?4 (logdiat = 0.81), respec- 
tively. From this we can derive that the RR Lyr surface 
distribution is flattened in the observed outer part, with 
b/a = diong/rfiat ^ 0.75. Complete data from the near- 
IR VVV survey will allow to study the distribution in 
the very central part of the bulge and to determine final 
parameters for the bulge RR Lyr population. 

6. CONCLUSIONS 

We have analyzed the data on 16,836 RR Lyr stars ob- 
served toward the Galactic bulge during the third phase 



located within 6° wide longitudinal (in red) and latitudinal slices 
(in blue). All points with logo! < 0.477 (in magenta) are common. 
There is no major difference between the three distributions, which 
further confirms the metal-uniform nature of the bulge RR, Lyr 
population. The steeper slope for the longitudinal slice indicates 
that the population is flattened along Galactic plane. 



of the OGLE project (years 2001-2009). After eliminat- 
ing low-amplitude stars, objects likely foreground and 
background sources, and GC members the sample in- 
cludes 10,472 RRab, 4608 RRc, and 78 RRd bulge vari- 
ables with derived mean V- and /-band brightness. For 
an additional 627 RRab, 151 RRc, and 2 RRd stars, all 
located close to the Galactic plane, we possess mean mag- 
nitudes in the / filter alone. 

Using color information for this type of standard can- 
dles, we find that the ratio of total to selective extinction 
Ri = Ai/E(y — I) is independent of color and equals 
1.080 ± 0.007. This result confirms the anomalous na- 
ture of the extinction toward the bulge. 

Our analysis clearly demonstrates that the bulge RR 
Lyr stars constitute a very uniform population, different 
to the majority of the stars in the central regions of the 
Milky Way. The photometrically derived metallicity dis- 
tribution is sharply peaked at f Fe/H],^^ = ~l -02 ± 0.18 
with a dispersion of 0.25 dex on IJurcsiM ()1995D metallic- 
ity scale. Based on the distribution of dereddened bright- 
ness we show that the RR Lyr population is elongated in 
the inner part, where it tends to follow the barred dis- 
tribution of the bulge RCGs. The elongated structure 
can also be noticed as a flat peak in the distance distri- 
bution. A possible explanation of this structure along 
the bar is that less populated old RR Lyr stars may feel 
gravitational forces from young stars in the Milky Way 
bulge. Farther off the Galactic plane (for |6| > 4°) the 
shape of RR Lyr population becomes spherical. We note 
that, due to small number statistics, we can not conclude 
the presence of a split in the population of bulge RR Lyr 
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stars, as is observed in RCGs. 

The distance to the Milky Way center inferred from the 
OGLE-III bulge RRab variables is Rq = 8.54 ± 0.42 kpc. 
This result is in agreement with Rn = 8 . 1 ± .6 kpc ob- 
tained from the OGLE-II data (Maiacss 2010) an d with 
a valu e oi Rq = 8. 15 ±0.49 kpc combined bv Genzel et al.l 
(|2010f ) from 11 different methods. It seems that the RR 
Lyr method suffers from uncertainties in the metallic- 
ity determination, and hence absolute brightness of the 
stars. 

In our work, we have also studied the density distribu- 
tion of the bulge RR Lyr stars as a function of distance 
from the Galactic center. We found a break in the dis- 
tribution at a distance of ~ 0.5 kpc indicating flattening 
toward the center. Taking into account completeness of 
the OGLE-III data, we estimate the number of RRab 
stars which should be detected within the bulge area of 
the near-IR VVV survey as (4-7) x lO'' stars. The un- 
certainty partially results from unknown density distri- 
bution within 0.2 kpc from the center. 

On going and future wide-field photom etric surveys, 
such as OGLE-IV, VVV, Pan-STARRS (jKaiser ct al. 
2002 1. and Large Synoptic Survey Telescope (Ivczic ct al. 
20081) should bring definitive answers to questions per- 
taining to the structure and evolution of all Galactic RR 
Lyr populations. For example, one of unsolved issues is 
whether the bulge and halo RR Lyr va riables form the 
same population. Mi nniti et al.l (|1998D suggested that 
the bulge RR Lyr stars could represent the inner ex- 
tension of th e Galactic h alo. The bulge pulsators do 
not show the lOosterhoB ()1939[ ) dichotomy observed in 
the halo variables. From the period distribution and 



period -amplitude diagram presented in iSoszvhski et al.l 
(|201lD it seems that probably all bulge RR Lyr stars 
belong to the Oosterhoff type I (Ool). Our work shows 
that the bulge population is very metal uniform and was 
likely formed on a short time scale. The Galactic halo 
was probably formed by at l east two distinct accretion 
processes (Mi ceh et al.ri2008[) . About 25% of the halo 
RR Lyr variables observed at distances 3-30 kpc ap- 
pear to be Oosterhoff type II (OoII) stars with a sig- 
nificantly different radial density profile to the remain- 
ing 75% of Ool objects. The Ool halo population has a 
power-law exponent of — 2.26±0.1 while the Poll compo- 
nent has a steeper slope with —2.88 ± 0.04 (jMiceli et al.l 
I2008D . The Ool stars are on average slightly more 
metal-rich (([Fe/Hjjgg) ~ —1.7 dex) than th e OoII stars 
(([Fe/H]jg5) 2.0 dex; ISesar et al.l [20Toh . Compari- 
son of the above numbers for the halo RR Lyr variables 
with those for the bulge objects obtained in this work 
does not allow us to answer the question on the common 
bulge/halo population. More data, in particular on RR 
Lyr stars residing in the central and outer regions of our 
Galaxy, are needed. 
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